function A = getA(E,N,alpha,gamma)
A = (E.^((gamma*(1-alpha)+1)/gamma)./(alpha*N.^(1/gamma))).^(1/(1-alpha));